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Abstract 

We consider the shear viscosity of a system of quarks and its ratio to the entropy density above 
the critical temperature for deconfinement. Both quantities are derived and computed for different 
modehng of the quark self-energy, also allowing for a temperature dependence of the effective mass 
and width. The behaviour of the viscosity and the entropy density is argued in terms of the 
strength of the coupling and of the main characteristics of the quark self-energy. A comparison 
with existing results is also discussed. 
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I. INTRODUCTION 



We address the problem of the viscosity of a Fermi system, in particular of quark matter. 
Much interest on this physics has been stimulated by the experiments carried out at the 
Relativistic Heavy Ion Collider (RHIC): these seem to indicate that above the critical tem- 
perature Tn the quark-duon plasma (QGP) behaves as a fluid with a quite small viscosity 

nnn 

(i.e. as an almost perfect fluid). This conclusion Ijj, y, yl, mostly based on the findings on 

nnriQ 

the f 2 coefficient measured 15|, la, 17|] in the multipole analysis of the angular distribution of 
the hadrons produced in an ultrarelativistic ion-ion collision, contrasts with a description of 
the QGP as a fluid of almost independent Landau quasi-particles which emerges from lattice 
calculations [8] though more recent findings substantially lowered also the lattice predic- 
tions. Also the substantial collective flow (the "elliptic" flow 



observed in these collisions 

appears to imply that the viscosity can not be that large |10l. I 111 . 12, [isl, [l^. In this work we 
shall examine this challenging question, conflning ourselves, at present, to the quark degrees 
of freedom only : hence our results will not be directly comparable with the indications of 
the RHIC data, but should be viewed as flrst a step in that direction, in the absence of the 
gluonic degrees of freedom. 

We will consider, beyond the shear viscosity rj, also the ratio rj/s, s being the entropy 
density of the system. Indeed rj/s, which for vanishing chemical potential basically governs 
the damping rate of sound waves propagating into the system [llj], has much physical signifl- 
cance: its temperature behaviour should determine the value Tc at which a phase transition 



occurs in the system [3|, [iSj. 



It should also be added that certain special supersymmetric = 4 field theories (dual to 
black branes in higher space-time dimensions) predict a lower limit for the viscosity /entropy 



density ratio, namely rj/s ^ I/Att (in units h = ks = c = 1) |16 |. 

After introducing, in Section [Tll the basic formalism to derive the shear- viscosity and the 
entropy density starting from a dressed quark propagator (and the corresponding spectral 
function) we present in Section IIIII our results and concluding remarks. 
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II. FORMALISM 

We base our investigation, along with many authors Q, 0, Q, Q, |20| . on the hnear 
response theory which, as far as the viscosity is concerned, leads to the Kubo formula : 



CO 

where the retarded Green's function at zero momentum is defined as 
n^M = -t I dte''^' / d^r < [T^,y{r,t),T^y{0,0)] > 



(1) 



(2) 



and the integrand is the correlator of the off-diagonal (x, y) element of the energy-momentum 
tensor T^jy at different space-time points. In Eq. ([2]) the brackets imply the thermal average 
of the commutator. 

From ([1]) and ([2]) the static {u = 0) viscosity follows: 



rj = r](uj = 0) = - —lmU^(uj] 



(3) 



=0+ 



We conveniently get from the Matsubara formalism by analytic continuation, accord- 
ing to the prescription {6 0) 



with 



(4) 



in the above Un = 27Tn//3, (3 = 1/T and is the r-ordered product. 

Concerning the structure of the canonical energy-momentum tensor we observe that for 
all the Lagrangians not displaying a derivative coupling among the fields (a typical example 
being the Nambu-Jona-Lasinio one) it turns out that 



Txy = - (V'72<9i^/' - <9i?/'72V') 



(5) 



ijj being the quark 
Following Ref. 



he 



d. 



17| . in the evaluation of Il{iuJn) one can stick to the first order of a 



ring diagram expansion, since higher orders (e.g. for a scalar or pseudoscalar interaction 



3 



coupling) vanish, in the chiral hmit ruq = 0, due to the trace of odd numbers of 7 matrices^. 
Hence, by inserting ([5]) into (jlj), one easily gets: 

1 y f d^p 

n(«^„) = -^2-/ J ^^^Px-Tr [725'(iwi + iw„,p)72S'(2aJi,p)] (6) 

where the trace is taken over spin, flavor and colour. 

Next, by introducing the spectral representation of the quark propagator, 

5(-..P)=r^^. m 

y_oo 27r luji - e 



with uji = {21 + l)n/l3, it is possible [17| to carry out the summation over the Matsubara fre- 
quencies via standard contour integral technique. One ends up with the following expression 
for the shear viscosity: 

where /(e) = l/(e^^^~'^^ + 1) is the thermal distribution for fermions; we notice that the 
viscosity 77 gets its major contribution from the surface of the Fermi distribution. 

In equation ([8]) the quark propagators are meant to be "fully" dressed within the appro- 
priate mean-field approach, depending upon the model Lagrangian one is referring to. Here 
we shall adopt a merely phenomenological approach, with a suitable Ansatz for the quark 
self-energy. 

A. The scalar case 

We consider first the case of a Lorentz scalar self-energy of the form 

E(p) = [M(p)-zr(p)]l (9) 
to be inserted into the (massless) fermion propagator 



S = — ^ • (10) 



^ A different situation would be encountered by considering series of ladder diagrams [.18|; however one 
should keep in mind that exchange-like diagrams (interaction inside a fermionic loop) are reduced with 
respect to the direct ones by the degeneracy of states, including the number of colors Nc- 
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In the above, M and F should be viewed as phenomenological functions of the four- 
momentum p, which will be specified later on. 

The spectral function associated with the propagator (ITUl) is then easily found to be 



rsign(e) iif + Mf + T^) 
p(e, p) = -I [^A(e, p) - SR{e, p) = 2 ^ ^ ^ > )_ 

where Sa and Sr are the usual advanced and retarded propagators. ^ 

After inserting (llip into ([H]) and performing the relevant traces one gets the following 
expression for the viscosity: 

v = 'ni + vii (12) 

with 

" T J (27r)^ 7-00 2vr ^ M M ; ^^^2 _ m2 + 12)2 + 4M2r2]2 ^ ^ 

and 

= - ^ y (2^ y_ 2^ - ^ (p2-M2 + r2)2 + 4M2r2 • (^3^) 

In the above = — p^, A^^c and Nf are the colour and flavour numbers. The temperature 
(and chemical potential) dependence is embedded into the Fermi distribution /(e). In the 
similar approach by Iwasaki et al. \y\ M and F were kept as positive, constant parameters; 
however with this choice rjjj, although generally smaller than rjj if integrated up to a cutoff 
momentum, is divergent. 

In order to ensure the convergence of rjjj and guided by simplicity arguments, we shall 
consider a constant parameter M, while for the width of the quasi-particle we use: 

A being a constant parameter as well. One should notice that the choice of the two parame- 
ters entering into the quasi-particle self-energy is not completely arbitrary, since on the basis 



of general arguments the spectral function must obey the following "sum rule" [2l| : 

1 de 

4^spiny ^-[p(e,p)7o] = l, (15) 



Notice that pip differs from the Lorentzian shape one would obtain in the non-relativistic hmit. 
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which is satisfied when the fermion propagator obeys a dispersion relation. 

Let us now consider the entropy density within the same model (we stick in this paper 
to the fj, = case); according to the customary field theory formulation ^] it reads: 

. =i^(TlnZ) (16) 



+ e'^'^+ 1 + e^"^- 
where V is the normalization volume and 

ujI (p) = p2 + - ± 2iTM . (17) 

By separating real and imaginary parts of uj±{p) it can be explicitly shown that s is real; 
one gets (in the hypothesis T < M): 



s = 2N,Nf I { In 



^ Ep ie p + cos Op] + Op sin ap 



+ 2e-^p ^ ' ()18) 

g-sSp _|_ 2e~^p cos Op + 1 



, e"^^'' + 2e"-^'' cos Op + 1 

(2vr) 

with Ep = PRe uj+{p), 5p = /?Im uj+{p). 
B. General self-energy 

Let us now turn to a more general structure for the quark self-energy: 

S(p) = ao(p)7° + ai(p) y • p + ^2(^)1 (19) 
p being the quark four-momentum. The corresponding propagator and spectral function can 



be more conveniently expressed making use of the customary projector operators 23|, l25l |: 



A.(p) = i(l±7o^^) (20) 



2 V 

where Ep = a/p^ + m'^ and m is the bare fermion mass. 
It is then easily shown that: 

joS-\p) = 70 - m - S(p)) = A-\p) A+(p) + Az'ip) A-(P) , (21) 

with 

A±^(po,p)=PoT[i?p + S±(p)]. (22) 
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In the above the self-energy (fT9i) has been expressed as well in terms of the projection 
operators fl20|) . with the condition 02 (p) = mai(p): 



7oE(p) = A+(p) - S_(p) A_(p) , (23) 

where 

^±ip) = aiip)Ep±aoip). (24) 

The functions ai{p) {i = 0,1) should be obtained on the basis of some microscopic calculation, 
thus entailing the modeling of the quark propagator in a suitable description. 

Attention should be payed to the definition of the retarded and advanced propagator, for 
which analogous definitions hold, e.g. with 

<^'^(^0,P)= ^.^'.(.), (25) 
Po T [Ep + '{p}] 

By definition the "retarded" self-energy must satisfy the condition: 

ImS^ < and ImS^ > (26) 

and 

^iip) = [^^{P)Y ■ (27) 
The spectral function can now be expressed as 

piPo, P)7o = P+{po, p) A+(p) + P-{po, p) A_(p) (28) 

where 

P±iPo, P = T ^ o ■ 29) 

[poT(i?p + ReS^(p))]' + [ImS^(p)f 

With these ingredients we can now turn to the evaluation of the shear viscosity and of 
the entropy density; the former is obtained from formula ([8]) and again can be split into two 
terms, 777 and rju, which read: 



NcNf f d^p r de pIpI . , . , ..2 i^r, ^ 

= ^ y ^ (1 - m) m [p,(e,p) - p_(e,p)] (30a) 
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and 



Vii 



NcNf 
T 



(27r)3 



-(l-/(e)) f{e)2plp^{e,p)p^{e,p) 
-oo 27r 



(30b) 



For what concerns the entropy density, we start from the general formula 



25| 



-2N,Nf 



(27r)3 



de dfie) 



Trspin {Im [In (-70^;^')] + Im (70^^) Re (^,j7o)}(91) 



,(27r) dT 

the trace being taken over the spin degrees of freedom. Having expressed all relevant quan- 
tities in terms of the projector operators, the trace of the logarithm can now be carried out, 
yielding 



s = -AN^Nf 
+arctg 



d^p 



de dfie) 



(2vr)3 (27r) dT 
ImS^ 



arctg 



-ImS^ 



e-{Ep + ReS^ 



e + Ep + ReS 



R 



^ -n9[e- {Ep + ReE^)] -ne[e + Ep + ReE^ + 



-ip+(e,p)[e - {Ep + ReS^)] - ^P-{e,p)[e + Ep + ReS^]| . 



(32) 



We now apply the above formalism to a schematic model for the quark self-energy. We 
stick to the chiral limit (m = 0) and start from the HTL self-energy 2J, |25| {p = |p|): 

p 



Sr^(e,p) = — (1 



e=FP, e + P 
m 



(33) 



2p ^ — V / 

then, by taking into account that at yU = the integrand in the shear viscosity is peaked at 
zero frequency, we consider the e — >^ 0+ limit (namely e — > i5, 5 being a positive infinitesimal) 
of the expression fl33l) : 



^HTL 



(0+,p) = T^vr— . 



(34) 



p 2p 

In the above formulas the temperature dependent mass parameter is M = AT, where (in 
the original HTL formula) = Cjg'^/S = g'^/Q for N^. = 3 colors, Cj being the Casimir 
invariant of SU{Nc) and g the strong coupling constant. 

The shear viscosity is then obtained by inserting into eqs. fISOap and (ISObp the following 
spectral densities: 

ttM^ 1 



P±(e,p) 



p 



M2 



+ 



2p 



2 • 



(35) 



In order to evaluate the entropy density fl32|) one can conveniently employ the identity 

dfie) _ dcTfie) 



dT 



de 



(36) 
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with 



af{e) = - [/(e) log /(e) + (1 - /(e))log(l - /(e))] ; (37) 



it allows to integrate by parts over the energy in (132|) . yielding the simple expression^ 



The function 7r^|p^(e,p) +p?.(e,p)} can be interpreted as the spectral density 26] 
for the entropy. Notice that in getting (!38|) the energy independence of the self-energy 
f lM|) was crucial; moreover we removed the contributions stemming from the derivative of 
the theta-functions in Eq. (!32l) . since they would have produced extra delta-peaks in the 
spectral density, leading to a double counting of the degrees of freedom (and altering the 
normalization of the spectral function itself). 

III. RESULTS AND CONCLUSIONS 

Using the formalism of the previous Section, we have evaluated the shear viscosity, the 
entropy density and their ratio with three different choices for the quark self-energy. For the 
simplest scalar self-energy [Eq. (Q] we first consider a constant M and F given by Eq. (IT^ 
with A = 0.5M; the choice for the latter was constrained by the sum rule for the spectral 
density, which, for A < 0.5M is only mildly violated (e.g. with the above value it differs 
from 1 by less than about 10% in the interval < p < 150 MeV). 

In Fig. [U we display the result obtained for rj/s with M = 300 MeV and 350 MeV, as a 
function of temperature for T > 140 MeV: the ratio grows like aT^, the coefficient a becoming 
smaller as M increases. We expect this model to be meaningful only in the regime of high 
temperatures, where it represents a gas of quasi-particles. It might be interesting to compare 
these results with the corresponding ratio obtained in a pion gas [^^l, also displayed in the 
figure. In this approach the shear viscosity is derived from a Boltzmann-Uehling-Uhlenbeck 
transport equation via the Chapman-Enskog method 28| to first order, while the entropy 
density is the local equilibrium bosonic entropy density. Since the differential cross section 
used in the numerical computation is the experimental one, these values of r]/s at zero 



The finite term vanishes since lim^^ioo ^^fi^) = 0. 
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FIG. 1 The ratio rj/s in the scalar model with M = 300 MeV (dashed Une) and M = 350 MeV 
(dot-dashed line) as a function of temperature for T > 140 MeV; the points at low temperatures refer to 



the pion gas (Courtesy of D. Davesne 



chemical potential are essentially model-independent. The two regimes (composite hadrons 
at low temperature, gas of quasi-particle quark states in the high T bath) qualitatively 
reproduce the analogous situation illustrated in Ref. [sl, where the quark-gluon phase is 
described within perturbative QCD. 

Within the same model, but taking into account a temperature dependence of the quasi- 
particle mass and width, the results for the ratio rj/s drastically change: indeed by heuris- 
tically assuming the conventional thermal free mass, M = 2ttT and, in eq.l fT^ . A = 6M, 
we obtained a T-independent rj/s ratio, which displays a 1/6^ behaviour, as it is illustrated 
in Fig. [2J We notice that for the largest value of h reported in the figure, our rj/s ratio is 
close to the AdS/CFT limit; moreover for h > 0.25 the shear viscosity becomes negative, 
an outcome which reflects the violation of the unitarity condition imposed by the sum rule 
( ITSl) : the latter indeed "protects" the positivity of r] and the balance between ReS and ImS: 
both grow linearly with T, but their ratio must be kept within well defined limits (with 
A = 0.25M the violation of the sum rule is smaller that 5%). 

Finally we considered the self-energy of Eq. (134|) , which exactly satisfies the sum rule (fT5|) 
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for any choice of M. As in the previous case, we found that both the shear viscosity and the 
entropy density raise proportionally to T^, so that their ratio is temperature independent. 
This can be seen by dimensional analysis since, sticking to fi = and with a fixed strong 
coupling constant (namely a T- independent g), the energy scale is set by the temperature 
itself. 

In Fig. [3] we report the quantity rj/T^ from Eqs. (JSOaj), (I30bl) and s/T^ from Eq. (l38ll 
as a function of A, while Fig. H] shows the corresponding ratio r]/s: the latter is compared 
with the same quantity obtained in the scalar case, with a temperature dependent mass, 
the b and A parameters being linked (at T = 0) by the relation b = \/\^87c. We find that 
in both cases r]/s decreases as 1/A^ and (since the entropy density is slightly varying with 
A) the dominant trend is due to rj. This l/F behaviour (F being the particle width) is a 
general feature of the microscopic calculation of transport coefficients, which stems from the 
analytic properties of the fermion spectral function 18|; likewise in kinetic theory the same 
behaviour is expected since l/F is related to the mean free path of particles. 

Also quantitatively the results we find are compatible with previous ones. For example 
at A = 0.65 {g ~ 1.6) we find rj/T^ = 6.7 to be compared with rj/T'^ ~ 7 found by [29^ for 



11 



f7/T'"3, s/T^3 




a quark plasma near the critical temperature. For the same A we get rj/s c:^ 0.8 which is 



four times larger than the one obtained in Ref. [13| for a purely gluonic system in a similar 
approach. Even by pushing our model to A = 1 we would find the somewhat larger value 
?7/s ^ 0.4: it seems that in the present approximation the quark contribution to rj/s is larger 
than the gluonic one (at least if we restrict ourselves to values of A for which the "HTL" 
approximation can be reasonably applied). 

In spite of the different Lorentz structure of the self-energy the two ratios displayed in 
Fig. m are quite similar. We can thus infer that the relevant feature which governs the ratio 
?7/s in the fermionic system is the temperature dependence of the effective mass and width 
of the quark, while the detailed (Lorentz) structure of the self-energy does not appear to be 
of much relevance. 

In summary, we have considered in this paper the rj/s ratio for a system of quarks with 
the purpose of investigating the impact on this quantity of the effective quark self-energy, 
which, although modeled on the basis of simplicity, still grasped important aspects of a 
realistic description. All calculations were carried out at zero chemical potential. Only for 
the case of a constant mass parameter we find that rj/s grows with the temperature, as 
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FIG. 4 The ratio jy/s versus A as obtained with the self-energy (p4)) (continuous Hne) and with the scalar, 
T-dependent model (dashed line); The AdS/CFT limit is also reported. 



expected for a limiting case of a dilute gas system [3|, ll5|] . 

On the contrary, by assuming mass and width parameters proportional to the temper- 
ature, no matter upon the details of the self-energy, the ratio rj/s becomes temperature 
independent, since both r] and s acquire a behaviour. 

Concerning the specific values obtained for this ratio, in our models it obviously depends 
on the coupling parameter (6 or A) adopted in the definition of the quark self-energy; we 
found a general behaviour ofrj/s, in agreement with previous findings. Our results are 
compatible with very small values of the ratio, but the fundamental requirement related to 
the sum rule obeyed by the quark spectral function sets some intrinsic limits on the value 
of the parameter itself and the AdS/CFT boundary is not reached. 

In the present paper we have shown that, within relatively simple models, one can ob- 
tain a realistic estimate for rj, although predictions are hindered by the arbitrariness of the 
employed coupling parameters. Actually, at the present stage no definite conclusion can be 
drawn about the shear viscosity of QGP: several microscopic models are yielding estimates 
of T]/s, but no calculation still exists which takes into account all the relevant degrees of 
freedom. Moreover further indications are needed both from experiment, beyond the infor- 
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mation provided by V2 [3Cf], as well as 
of this quantity in the hadronic phase 



ro m hy dro-calculations 
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31 



32] • Also the behaviour 



34l | deserves further investigation in order to get 



quantitative information on the deconfinement phase transition. 
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